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Abstract 

Bohr's model agreed with the hydrogen spectrum results, but did not agree with the spectrum 
of Helium. Here we show that Bohr's model-based methods can calculate the experimental value 
(-79.005 eV) of Helium ground state energy correctly, we suppose the orbital planes of the two 
electrons are perpendicular to each other. By a computational method, we calculate the Coulomb 
force among the particles, and the number of de Broglie's waves contained in the short segment 
at short time intervals. Our results demonstrate that two electrons of Helium are actually moving 
around, not as electron clouds. 
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In 1913, Niels Bohr postulates the Bohr's model which agreed with the observed hydrogen 
spectrum [l(]. In 1923 Louis de Broghe suggested that electrons might have wave aspect and 
its wavelength A is equal to h/mv, where h is Plank's constant (= 6.62606896 x 10"^'^ 
Js) and m is the electron mass (= 9.1093826 x 10^'^^ kg). In 1927, Davisson and Germer 
experimentally confirmed de Broglie's hypothesis [2]. Recently the results of the two-slit 
experiment of an electron showed its wavelike properties 

In the Bohr's theory, the circular orbital length is equal to a integer times the wavelength 
of the electron, so we have, 27rr = n x h/mv. In several phenomena, the Bohr's model 
p..ov,de. good accu..acy Q Q. Q. But th. „,ode. could not explain aboot the spin of the 
electron and the two-electron atoms such as the helium exactly. Because of such problems, 
the Bohr's theory was replaced by the quantum mechanical theory based on the Schroedinger 
equation in 1920's. The solution of the Schroedinger wave equation showed that the orbital 
angular momentum of the electron in the ground state of the hydrogen atom is zero. And 
the solution also showed the relation between L (the total angular momentum quantum 
number) and M (the z component of L) 

M = -L,-L + l,-L + 2, ... , L-2,L-1,L (1) 

Eq. (pP) and the results of the Stern-Gerlach experiment indicated that an electron has 
spin angular momentum. But we can not visualize the electronic motion concretely. The 
Coulomb potential may be infinitely negative, when the electron is close to the nucleus. And 
by equating the angular momentum of the spinnig sphere of^he electron to i^/i, the sphere 
speed leads to about one hundred times the speed of light [7] . To solve the above problem, 
we try going back to the Bohr's model. 

The helium atom has two electrons and a nucleus of charge -|-2e. The Schroedinger 
equation for the helium can not be precisely solved. Using the perturbation method, we 
cariget t he g round state energy of the helium, which is very close to the experimental value 



ll| . But the calculation of the high order correction is very difficult. 
In this paper, we try to calculate the ground state energy of the helium atom using the 
new theory based on the Bohr's model, and check if the calculation value is equal to the 
experimental value -79.005 eV. 

First, suppose we have one model as shown in Fig. [H In that model, two electrons of the 
helium are on the opposite sides of the nucleus and moving on the same circular orbital. 
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FIG. 1: One schematic model of Helium in which two electrons are moving on the opposite sides 
of the nucleus. 



Equating the centrifugal force to the Coulomb force, we have 



(2) 



r Airer'^ 47re(2r)2 

where r is the circular orbital radius (m), e is the electron charge (= 1.60217653 x lO^^^C), 
and e is the permittivity of vacuum (= 8.854187817 x lO^^^-^^a ) The circular orbital length 
is supposed to be an integer times the wavelength of the electron, we have 

h , . 

27rr = X n (3) 

mv 

The total energy E of the electrons is the sum of the kinetic energy and the Coulomb 
potential energy, so 

E = 2x 2 X \ — - (4) 

Solving the above three Eqs. ([2lll]), the ground state energy (n=l) is -83.33 eV. The value is 
lower than the experimental value -79.005 eV. In this model, as the two electrons are on the 
same orbital, it is possible that de Broglie's waves of the two electrons interfere with each 
other and their motions are affected, so Eqs. (ISlHj) may not be satisfied. Whether the wave 
is the electron itself or by the field change around the electron can not be ascertained. 

To avoid such problems, we suppose another model as shown in Fig. [2] and [3l In that 
model, the planes of the two same-shaped orbitals are perpendicular to each other. As any 
point on the electron 1 orbital is at the same distance from the points on the both-side (± 
z) electron 2 orbital, de Bloglie's wave of the electron 2 may interfere with itself and vanish 
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FIG. 2: Schematic model of Helium. Two same-shaped orbital planes are perpendicular to each 
other. This figure shows one quarter of the orbitals. Electron 1 starts at (rl, 0, 0), while electron 
2 starts at (-rl, 0, 0). 
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FIG. 3: Schematic model of Helium. Two electrons have moved one quarter of their orbitals. 
Electron 1 is crossing y axis perpendicularly, while electron 2 is crossing z axis. 

on the electron 1 orbital. So the motion of electron 1 may not be affected by the wave of 
electron 2. 

Here we investigate how the electrons of the helium are moving by calculating the 
Coulomb force among the two electrons and the nucleus at short time intervals. The com- 
puter program (class filename: MathMethod) written in the JAVA language (version 1.5.0) 
to calculate the electronic orbital of the helium is shown in Supplementary Methods. 

As shown in Fig. [2] and [3], the helium nucleus is at the origin, the orbital plane of the 
electron 1 is the xy-plane, and that of the electron 2 is xz-plane. The electron 1 initially 
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at (rl, 0, 0) (Fig. [2]) moves one quarter of its orbital to (0, r2, 0) (Fig. [3]), while the 
electron 2 initially at (-rl, 0, 0) moves to (0, 0, r2). Meter and second are rather large 
units for measurement of atomic behavior, here we use new convenient units MM (1 MM 
= 1 X 10-1^ meter), SS {1 SS = 1 x IQ-^^ second) and MM/SS (1 MM/SS = 1 x 10^ 
meter /second). 

In this program, we first input the initial x-coordinate rl (in MM) of the electron 1, and 
the absolute value of the total energy E (in eV) of the helium. From the inputted value, 
we calculate the initial velocity of the electron. And at intervals of 1 SS we calculate the 
Coulomb force among the two electrons and the nucleus. When the electron 1 is at (xx, yy, 
0), the electron 2 is at (-xx, 0, yy) (in MM). Change MM to meter as follows; x (m) = 
XX xlO"^"^. y (m) = yy xlO~^^. So the x component of the acceleration (m/sec^) of the 
electron 1 is 

2e^x X (2a;) 
T + \ \. ..3 (5) 



47re(rm) (x^ + y'^) 2 A7re{rm) {Ax"^ + 2y'^) 2 
where the first term is by the Coulomb force between the nucleus and the electron 1, and 

the second term is by the force between the two electrons. Considering the helium nuclear 

mass, we use here the reduced mass (rm = 9.1069044 xlO~^^ kg) except when the center of 

mass is at the origin. 

In the same way, the y component of the acceleration [m/sec^) is 

26^-0 X y , . 

+ (6) 



47re(rm) (x^ + y"^) 2 47re(rm) (4a;^ + 2y'^) 2 
Change m/sec^ to MM/SS"^ using the next relation 

Im/sec^ = 1 X IQ-^^MM/SS"^ (7) 

Based on that calculation value we change the velocity vector and the position of the elec- 
trons. We suppose electron 1 moves only on the xy-plane, so the z component of the 
acceleration of the electron 1 is not considered. We also calculate de Broglie's wavelength 
of the electron from the velocity (A = h/mv) at intervals of 1 SS. The number of that wave 
(A in length) contained in that short movement section (The sum of them is WN) is 

VVX^ + VY^ X 10-^^ , , 

h (-8) 

(rm 

where {VX, VY) are the velocity of the electron 1 (in MM/ SS), the numerator is the move- 
ment distance (in meter) for 1 SS. the denominator is de Broglie's wavelength (in meter). 



TABLE I: Results of rl and WN (Number of de Broglie's waves) in which y component of electron 
1 velocity in Fig. [3] is zero at various energy levels of Helium. WN x 4 is the total number of de 
Broglie's waves contained in one round of the orbital. 
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Here, the estimated electronic orbital is divided into more than one hundred thousand for the 
calculation. When the electron 1 has moved one quarter of its orbital and its x-coordinate 
is zero (Fig. [3]), this program checked the y-component of the electron 1 velocity (last VY). 
When the last VY is zero, two electrons are periodically moving around the nucleus on the 
same orbitals as shown in Fig. [2] and O So, only when -0.0001 < VY < 0.0001 (MM/SS) is 
satisfied, the program displays the following values on the screen, rl, VY, preVY (VY ISS 
ago), and (mid)WN (the total number of de Broglie's waves contained in one quarter of the 
orbital). The initial inputted x-coordinate of the electron 1 is automatically increased by 1 
MM per above calculation to +100 MM. 

Table [T] shows the results in which the last VY is the closest to zero. Fig. H] graphs the 
results in Table I. As shown in Table I and Fig. HI when the total energy of the helium (E) 
is -79.005 eV (the experimental value), WN x 4 is just 1.00000. This results demonstrate 
that two electrons of the helium are actually moving around the nucleus on the orbitals of 
just one de Broglie's wavelength as shown in Fig. [2] and [31 
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FIG. 4: Plots of the number of de Broglie's waves contained in one orbital at various energy levels 
of Helium 



Experimentally, the spin magnetic moment (/i) is known to be 

g(3S 



h 



,/? = — 

2m 



(9) 



where g, the spin g-factor, is 2 5* is spin angular momentum, and /? is the Bohr 
magneton (= 9.27400949 x 10~^^^). If we suppose that the orbital angular momentum 
(5*2) of the electron in the ground state of the hydrogen atom is ±h (here, spin is zero) 
as in the Bohr's model and its g factor {g2) is 1, the magnetic moment of the electron 
fi2 = —{g2)[3{S2)/h is equal to of Eq. ([9]). In this state, the above-mentioned problems 
of infinitely negative Coulomb potential and spinning sphere speed do not occur. As the 
magnetic moment of the electron is not changed, the splitting pattern of the energy levels 
in the magnetic field (Paschen-Back Effect) is also not changed. 

In this study, the two main forces that affect the electronic motion are supposed to 
be the Coulomb force and the force by de Broglie's waves. By a simple calculation, the 
magnetic force between the two electrons is too small in comparison to the Coulomb force. 
In the standard model, as the helium atom has spin-up and spin-down electrons, it seems 
to generate no magnetic field. But to be precise, in all areas except in the part at just 
the same distance from the two electrons, magnetic fields are theoretically produced by the 
electrons even in the standard helium model. So as the electrons moves, they lose energy by 
emitting electromagnetic waves. The same thing can be said about this new helium model. 
In this model, the length of both orbitals crossing perpendicularly is just one de Broglie's 
wavelength. Probably the stability of de Broglie's waves are related to the stability of the 
electronic motions. 
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It is a very meaningful and astonishing thing that the calculation value by this new 
method based on the Bohr's model is just equal to the experimental value. As the Bohr's 
model applies to the hydrogenhke ion He+, the second ionization energy (54.418 eV) of 
the helium is correctly calculated using the Borh's theory. So, the first ionization energy 
(24.587 eV) can also be gotten by the equation 79.005 eV - 54.418 eV = 24.587 eV. This 
fact means the atomic electrons are actually moving around the nucleus, not existing as 
electronic clouds, it is possibile that the basic conception of the quantum mechanics such 
as the uncertainty principle will be greatly changed. 
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Supplementeiry Methods; JAVA computer program to calculate the electronic 
orbital. 

In this program, we first input the initial x-coordinate rl (in MM) of electron 1, and 
the absolute value of the total energy E (in eV) of Helium. From the inputted values, this 
program outputs the y component of electron 1 velocity in Fig. 3, and WN (the number of 
de Broglie's waves contained in one quarter of the orbital. 



import Java. util. Scanner; 

class MathMethod { 

public static void main (String[ ] args) { 

Scanner stdin = new Scanner ( System.in ) ; // input rl and |E | 

System.out.println ("rl between nucleus and electron 1 (MM)? "); 
double r = stdln.nextDouble (); 

System.out.println ("total energy |E |in the herium (eV) ? "); 
double E = stdln.nextDouble (); 

double me = 9.1093826e-31 ; double mp = 1.67262171e-27; 
double mn = 1.67492728e-27 ; double nucle = 2*mp + 2*mn; 

/ /rm = reduced mass of an electron 
double rm = (2*me*nucle) / (2*(2*me+nucle)); 

double pai = 3.141592653589793; double epsi = 8.85418781787346e-12; 

double h = 6.62606896e-34; double ele = 1.60217653c-19; 

for (int i = 1; i <100; i++ ) { // repeat until rl=initial rl + 100 

/ / calculation of initial VY from E and rl 
double poten = - (2*ele*ele*2) / (4*pai*epsi*r) + (ele*ele) / (4*pai*epsi*2*r); 

/ / vya = total E-potential energy 
double vya = - (E*1.60217646e-19) - poten*1.0el4; 
if (vya >0) { 

/ / vyb=velocity from kinetic energy 

double vyb = Math.sqrt(vya/me); 

double VY = vyb*1.0e-8; // change m/sec to MM/SS 

double prexx = r; double VX = 0.0; double WN = 0.0; double preyy=0.0; 
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double XX, yy, vk, preVY, preWN, midWN; 
do { 

XX — prexx + VX; yy = preyy + VY; / /electron 1 position after ISS 

preVY = VY; preWN = WN ; 

vk = VX * VX + VY * VY; // calculation of WN from VX,VY 

// WN = WN + Eq(8) 

WN = WN + (rm*vk*1.0e-6) / h; 

/ / calculation of VX,VY from Coulomb force 
double ra — Math.sqrt (prexx * prexx + preyy * preyy); 
double rb = Math.sqrt (4.0 * prexx * prexx + 2.0 * preyy * preyy); 

/ / change MM to meter 

ra = ra * l.Oe-14; rb = rb * l.Oe-14; 

prexx = prexx * l.Oe-14; preyy = preyy * l.Oe-14; 

double ac = (2 * ele * ele) / (4 * pai * epsi * rm); 

// VX = VX + Eq(5) * 10-=^° 
VX = VX + 1.0^30 * ac * prexx * (-1.0 / (ra*ra*ra) + 1.0 / (rb*rb*rb)); 

// VY = VY + Eq(6) * 10-30 
VY = VY + l.Oe-30 * ac * preyy * (-1.0 / (ra*ra*ra) + 0.5 / (rb*rb*rb)); 
prexx = xx; preyy = yy; 

} while ( XX >0 ); / /repeat above until electron 1 arive at y axis 

if (VY >-0.0001 && VY <0.0001) { // last VY condition 

System.out.print ("rl: " + r + " "); 
System.out.printf ("VX: % .5f ", VX); 
System.out.printf ("VY: % .5f ", VY); 
System.out.printf ("preVY: % .5f ", preVY); 

midWN = (preWN + WN) / 2; System.out.printf ("midWN: % .5f \n", midWN); 
} 

}r = r + 1; 
}}} 
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